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ABSTRACT 



Emission mechanisms of the shell-type supernova remnant (SNR) RX J1713.7- 
3946 are studied with multi-wavelength observational data from radio, X-ray, 
GeV 7-ray to TeV 7-ray band. A Markov Chain Monte Carlo method is employed 
to explore the high-dimensional model parameter space systematically. Three 
scenarios for the 7-ray emission are investigated: the leptonic, the hadronic and 
a hybrid one. Thermal emission from the background plasma is also included 
to constrain the gas density, assuming ionization equilibrium, and a la upper 
limit of about 0.03 cm -3 is obtained as far as thermal energies account for a 
significant fraction of the dissipated kinetic energy of the SNR shock. Although 
systematic errors dominate the x 2 of the spectral fit of all models, we find that 1) 
the leptonic model has the best constrained model parameters, whose values can 
be easily accommodated with a typical supernova, but gives relatively poor fit 
to the TeV 7-ray data; 2) The hybrid scenario has one more parameter than the 
leptonic one and improves the overall spectral fit significantly; 3) The hadronic 
one, which has three more parameters than the leptonic model, gives the best 
fit to the overall spectrum with relatively not-well-constrained model parameters 
and very hard spectra of accelerated particles. The uncertainties of the model 
parameters decrease significantly if the spectral indices of accelerated electrons 
and protons are the same. The hybrid and hadronic models also require an 
energy input into high-energy protons, which seems to be too high compared 
with typical values of a supernova explosion. Further investigations are required 
to reconcile these observations with SNR theories. 
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Subject headings: radiation mechanism: non-thermal — ISM: supernova rem- 
nants — cosmic rays — Gamma rays: ISM 



Introduction 



Supernova remnants (SNRs) are widely thought to be one important kind of cosmic 
ray (CR) sources in the Galaxy (lAharonian et al.ll2004l ). The most direct evidence comes 
from high-energy 7-ray emission from SNRs. Generally there are two types of scenarios for 
production of high energy 7-rays: the leptonic (via inverse Compton scattering of background 
photos by relativistic electrons) and hadronic (via decay of neutral pions produced by elastic 
collisions of relativistic ions with ions in the background plasma) origins. Understanding 
which of these two scenarios is dominant in specific sou rces is very im portant for the search 



of CR nuclei sources and the study of CR acceleration (jGabicil 12008 



Usually it is difficult to distinguish the leptonic model and hadronic model just with 
the high energy 7-ray data alone. Multi-wavelength observations of photon emission from 
SNRs can provide us key information about the radiation mechanism. Shell-type SNR 
RX J1713. 7-3946 is one of the most widely studied SNRs with perhaps the best m ulti- 
wavele ngth observations. The observational da ta span from radio (jLazendic et al.ll2004T). in 



2003 



2000 



Cassam-Chenai et al. 



Enomoto et al. 



2002 



frared (IBeniamin et al.ll2003l: lAcero et alil2009l), X-ray jKovama et al] 



ray emissio n detected by Suzaku (ITanaka et al 



1997 



Uchivama et al. 



2004h. GeV^-rav (Ubdo et al.lboilh , to TeV 7-ray band jMuraishi et al 



Aharonian et al. 12006). Re cent observations, especially the X- 



20081 ) and TeV 7-ray emission measured 



by HESS (jAharonian et al.l 120071 ). give the energy spectra and images of this SNR with 
yery high quality, wh i ch makes detailed modelings of the emission mechanism plausible 



(Morlino et a' 



Ellison et al 



2009a: 



2010 



Fan et al. 



Fang et al.l 120091 ; Fan et al.l l2010at IZirakashvili fc Aharonian 



2010b|). The newly reported data from Fermi flAbdo et al. 



2010 



20 111 ) also set strict constraints on the nature of the radiation from this SNR. 



Basic results of recent studies of this SNR may be summarized briefly as the follow- 
ing. The wide range TeV 7-ray spectrum favors a hadronic origin o f the high energy 



emission (lAharonian et al.l 120061 : iDrury et al.l l2009t iMorlino et al.l l2009al : iFang et al.l 12009 



Berezhko fc Volkil2010h . This scenario is also in line wi th the long standing view that SNRs 
are the most important CR accelerators (lAxfordlll98ll ). However, there is a strong correla- 
tion between the X-ra y image and TeV 7-ray image, favoring a le ptonic o rigin of the multi- 



wavelength emission ( lAharonian et al.l l2006t lAcero et al.l 120091 ) . iPlagal (120081 ) also claims 



that the lack of spatial correlation between 7-rays and the molecular cloud in the vicinity 
of SNR RX Jl 713. 7-3946 argues against the hadronic scenario. Furthermore the lack of 
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thermal line emission on th e X-ray spectrum sets an u pper limit on the ambient plasma 
density of about 0.02 cm -3 (jCassam-Chenai et al.l 12004 ) . which implies a very high energy 
content of accelerated protons from the supernova explosion. The hadronic model actually 
has a proton acceleration efficiency more than 4 orders of magnitude hig her than the elect ron 
acceleration efficiency, which corresponds to a rather extreme scenario ( jButt et al.l 120081 1 . It 
has been shown that the leptonic model can reproduce the multi -wavelength data and the 
model parameters can b e easily accommodated by typical SNRs (jZirakashvili fc Aharonian 
20101 ; Ellison et al.ll2010l ) though the overall fit to the data is relatively poor for the simplest 
cases. The spectral fit can be improved by considering; d etails of electron acceleration near 



the high energy cutoff (ILiu et al. 



20081 : iFan et al.ll2010al ]bl). In fact, the cutoff of the TeV 



spectrum at a few tens of TeV favors the leptonic model since production of photons at even 
higher energies through inverse Comptonization of the cosmic microwa ve radiation b y TeV 
electrons is in the Klein-Nishina regime and therefore very inefficient f lGabicil 120081 ). The 
observed decay of bright X-ray filaments with a width of ~ 0.1 lyr on a timescale of ~ 1 
year, on the other hand, can be attributed to fast diffusion of high energy electrons in a 
weak magnetic fi eld away from interm i ttently formed regions of relatively high accelerated 
electron density (jUchiyama et al.l 120071 ; ILiu et al.ll2008l ). 



It is evident that even with both theoretical and observational advances on this source 
recently, the nature of the TeV emission is still inconclusive. The multi-wavelength ob- 
servational data justifies a systematic modeling of the emission spectrum. Moreover, the 
thermal emission needs to be taken i nto accoun t to get more quantitative constraints. This 
paper focuses on these two aspects. Fan et al.l ( 12 01 Obi ) studied the goodness-of-fit for var- 
ious physically motivated leptonic models and found that the diffusive shock acceleration 
and the stochastic acceleration give comparably good fits. In this work we generalize this 
analysis by including the hadronic component. We employ the Markov Chain Monte Carlo 
(MCMC) method to constrain the model parameters, and investigate the full, correlated pa- 
rameter space systematically. The non-thermal spectra of CR electrons and/or protons are 
parameterized in the simplest way, i.e., a power law with a high-energy cutoff. The thermal 
bremmstrahlung radiation and line emission of the background plasma are also taken into 
account in the fit. We consider three scenarios for the 7-ray emission, the purely leptonic 
model, the hadronic model and a hybrid model where the number of model parameters is re- 
duced by requiring the spectral parameters of CR protons and electrons are identical except 
the normalization (see Sec. 2. below). The fitting results are presented in Section 2 and 
show consistency with previous studies. The conclusion is drawn in Section 3, where we also 
discuss possible future researches necessary to improve our understanding of this source. 
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2. Fitting results 



In this section, we use the MCMC technique to constrain the model parameters. The 
MCMC method is well suitable for high dimensional parameter space investigation. The 
Metropolis-Hastings algorithm is used when sampling the model parameters. The probabil- 
ity density distributions of the model parameters can also be simply approximated by the 
number density of the sample poi nts. A bri e f intro duction to the basic procedure of the 
MCMC sampling can be found in Fan et al. (feoiOb ). For more detai ls about the MCMC 
method, please refer to Neal ( 1993 ); Gamerman (jl997 ); Mackay ( 2003 ). 



We also discuss implications of model parameters from the best fits to multi-wavelength 
data of SNR RX J1713. 7-3946 for three scenarios of the 7-ray emission. In all these scenarios, 
the radio to X-ray emissions are generated through synchrotron of relativistic electrons. The 
high energy 7-rays are produced with different mechanisms. The basic physical parameters 
of SNR RX J1713.7 -3946 are adopted as: Age T life « 1600 yr, Distance d « 1 kpc, and 
Radius R ~ 10 pc f lWang et al.lll997l ). and we assume a uniform emission sphere with a 
radius of R in deriving related quantities. Alt hough the errors o f the Fermi data are large, 
we still include these data in the spectral fits ( lAbdo et al.ll201ll ). The procedure described 
in this paper can be applied to future observations with improved data to evaluate different 
emission models. 



2.1. Leptonic scenario 



In the leptonic scenario the 7-ray emission is produced through inverse Compton (IC) 
scattering of energetic electrons off the background radiation field, including the interstel- 
lar infrared, optical radiation, and the cosmic microwave background (CMB). The energy 
spectrum of accelerated electrons is prescribed as F e (E) oc E~ ae exp [— (E/E^) Se ^ , where E, 
a e , El are the electron energy, power law spectral index, high-energy cutoff energy, respec- 
tively, and S e describes the sharpness of this cutoff. The normalization is given through the 
total energy of electrons above 1 GeV, W e . The synchrotron radiation also depends on the 
magnetic field strength B. 

The interstellar radiation field (ISRF) other than the CMB may be important for the 
calculation of IC 7- ray spectrum. The inclusion of ISRF has been proposed to improve the fit 
to the HESS data (IPorter et al.ll2006l ). However, given the new data of X-ray by Suzaku and 
TeV 7-ray by HESS, it was shown that only if the intensity of ISRF is artificially boosted 
by more than one order of magnitude, the goodness-of-fit can be improved significantly 
( ITanaka et al.ll2008l ; iMorlino et al.ll2009a| ). In this work the ISRF is adopted as that given 
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by lPorter et al.l (120061 ) at a distance of 7.5 kpc from the Galactic center and in the equatorial 
plane of the Galactic disk (IMoskalenko et al.ll2006l ). Our results are not sensitive to details 
of the ISRF. 

Thermal X-ray emission mostly depends on the density of the shocked interstellar 
medium (ISM) nisM and the temperature of background electrons T e . Depending on the 
effect of accelerated particles on the shock structure , the density is a factor of a few higher 
than the density in the un-shocked upstream region (IBerezhko fe Ellisonlll999l ; I Warren et al. 
20051 ). The e lectron tempe r ature due to Coulomb c ollisional energy exchange with ions is 
estimated by lliughes et all feoooh : iFan et all (boiOah 



T P > 2.1 x 10' 



life 



1600 yr 



2/5 



™ISM 



cm 



2/5 



Ti 



1.3 x 10 8 K 



2/5 



where Tj is the temperature of background ion s. Tj is estimated to be higher than ~ 1.3 x 10 8 
K if the background is heated by the shock (|Fan et al.l l2010ai ) . For rijsM ~ 0.02 cm -3 , T e 
should be higher than 4 x 10 6 K. iDrury et al.l (120091 ) alternatively proposed that the post- 
shock region temperature could be reduced significantly in the case of large Mach number 
of the shock and effective particle acceleration. If this is the case, as we will show below, 
the constraint on the density and therefore the relativistic proton energy in the hadronic 
and hybrid models will be les s strict. However, there are still signifi cant uncertainties in the 
plasma heating downstream (jGhavamian et al.l 120071 : iGabicil 120081 ) , and it was also argued 
that such an extreme condition was not e asy to meet and the post-shock plasma should be 
heated more strongly (lEUison et al.ll2010l ). Since there is no direct constraint on T e due to 
the lack of thermal emission, instead of including T e in the MCMC fit, we take some typical 
values of T e . For the sake of simplicity, we assume that the background electrons are in 
ionization equilibrium with the background ions and use the Raymond- Smith plasma code 
to ca lculate the thermal emission for given nisM, T e and the metalicity ([Raymond fc Smith 
19771 ). The emission includes the recombination, bremmstrahlun g, two photo n process, and 
line emissions. The chemical abundance of the ISM is taken from lAllenl (I1973I ). For emission 
lines lying in the energy range of X-ray data (0.5 — 33 keV), we convolve the model spectrum 
with a Gaussian ener gy spread function, w hose width is adopted as the characteristic energy 
resolution of Suzaku (IKoyama et al.l 120071 ) . 



In total there are 6 free parameters, a e , W e , 5 e , B, and nisM in the lepton model. 
The 1-dimensional (1-D) probability distributions and 2-dimensional (2-D) confidence regions 
(at la and 2a confidence levels) of the model parameters, and the best-fit spectral energy 
distribution (SED) of the source are shown in Fig. [TJ The best fit model parameters 
correspond to the peak of the 1-D probability distributions. The spectral parameters in 
the leptonic scenario are well constrained except for nisM, whose 2a upper limit is well 
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Fig. 1. — Left: 1-D probability distribution of the parameters in the leptonic model; 
Middle: 2-D confidence contours of the parameters. The contours are for 1 and 2 a lev- 



els; R ight: th e best fit to the sp ectral energy distribution (SEP ) from radio ( lAcero et al. 



2009h X-rav JTanaka et "all 120081 ). GeV 7-ray (Ubdo et alJboilh and TeV 7-ray observa- 
tions ( Aharonian et al. 2007h . The background electron temperature is 10 7 K. The thermal 
emission component indicated by the dotted line corresponds to the 2a upper limit of 0.007 
3 for m SM . 



cm 



determined. In this calculation we set T e = 10 7 K. The only parameter sensitive to T e is 
nisM in the leptonic scenario. The 2a upper limit of nisM is 0.007 cm -3 for T e = 10 7 K. The 
dotted line in the right panel of Fig. [1] indicates the thermal emission for these parameters. 
Since the ISM density will be more essential for the discussion of the hybrid and hadronic 
models, we will discuss the T e dependence of these results in detail in the subsection 2.4. 



For the 2-D confidence regions of the parameters, we only show combinations with 
relatively large correlation. There are very weak correlations among a e , E^, W e and B. 
The weak correlation between a e and W e is mostly due to the facts that electrons near the 
high-energy cutoff are well-constrained by observations and low-energy electrons contribute 
the most to W e for a e > 2. The correlation between 5 e and E\ is caused by the well observed 
spectral shape in hard X-rays and TeV 7-rays. The combination of the X-ray and 7-ray data 
helps to determine the model parameters. In the following we will see that X-ray data alone 
lead to poorly constrained and highly correlated parameters. This is an example to show 
the importance of global fit to the multi-wavelength data. 

The parameters and \ 2 values of the best fit model are compiled in Tables [1] and |2l 
respectively. Since 7Zism is not well constrained, its 2a upper limit instead of the best fit 
value is listed in Tabl e U The best-fit parameters are c onsistent with previous studies 



( 1 Aharonian et al.l 120061 ; iLiu et al. 



2008 



Tanaka et al.l 120081 ) and can be readily accommo- 



dated with typical SNRs (lEllison et al.ll2010l ). The overall \ 2 °f the fit is relatively large 
with the reduced x 2 ~ 466.9/232 = 2.01. Such a high value of \ 2 shows that systematic 
errors dominate. This is not surprising given the relatively complex structure of the SNR, 
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uncertainties related to the particle acceleration process, and our rather simple prescription 
of the emission model. The systematic errors actually dominate in all emission models, which 
is typical for modeling of astrophysical observations of relatively complex phenomena. The 
X-ray and TeV 7-ray data contribute the most to the overall x 2 - Especially for the TeV 7-ray 
data, the x 2 value is 149 for 27 data points, corresponding to an average residuals about 
2.3<r. That is to say this simple leptonic model actu ally can not fit the HESS data well. 



This is a well-known resul t in previous stud ies ( e.g.. lAharonian et al.l 120061 ; iTanaka et al. 



2008 



Morlino et al.l l2009al ; iFang et al.l 120091 ) . In iLiu et al.l (120081 ) the authors proposed a 
stochastic acceleration model to generate the electron spectrum with sub-exponential cutoff 
(8 e = 0.5) to better fit the HESS data. However, in such a case the fit to X-ray data be- 
comes worse. The X-ray data actually favors super-exponential cutoff instead (with 5 e = 1.2 
in this p urely leptonic fit). The fit may be improved in some detailed leptonic models, as 
shown in iFan et al.l ( l2010bl ). though systematic errors still dominate. 



Table 1: Fitting parameters for T e 
correspond to 2a confidence level. 



10 7 K. Errors are la statistical uncertainties; limits 
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Table 2: Best-fit x 2 values for each set of data and the total reduced one, for T e = 10 7 K. 
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X-ray 


GeV 


TeV 


reduced 


leptonic 


2.07 


312.0 


3.4 


149.4 


466.9/232 


hadronic 


0.13 


291.7 


1.9 


40.6 


334.4/229 
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0.10 
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1.9 


40.5 


334.4/230 


hybrid 


0.02 


305.4 


20.0 


109.3 


434.7/231 



2.2. Hadronic scenario 



In this subsection we discuss the model with a predominantly hadronic origin of the 7- 
rays. The spectrum of the accelerated protons is assumed to be F P (E) oc E~ ap exp[— (E/E^) Sp ] 
with 8 P = 1, which gives acceptable fit to the TeV data. The normalization is fixed using 
the total kinetic energy of protons with th e energy E > 1 GeV . For the hadronic 7-ray pro- 
duction we adopt the parameterization of iKamae et al.l ( 120061 ). With the additional three 
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parameters, a p , E?, W p , we have 9 parameters in total. Considering the synchrotron cooling^ 
of high energy electrons, we also introduce a spectral break to the overall electron distribu- 
tion.]. The break energy, at which the synchrotron cooling time is equal to the lifetime of the 



remn ant, is determined with E hr « 7.8 x 10 6 ( J B//iG)- 2 (T life /1600 yr)~ 1 GeV flTanaka et al. 



20081 ). For E e < E^ r the power- law index is a e , and for E e > E^ the power- law index is 



a e + l. Since 2]if e is taken as 1600 years, E\> r is not a free parameter. 

The 1-D probability distributions and 2-D confidence contours of the model parameters, 
and the SED of the best-fit model are shown in Fig. [2J The parameters and \ 2 values of 
the best fit are listed in Tables [T] and |5J respectively. Still we adopt T e = 10 7 K in this 
calculation. Compared with the leptonic m odel, a much stronger m agnetic field is inferred, 



which is consistent with previous studies (IBerezhko fc Volki 120061 ) and the interpretation 



of the obs erved X-ray surface b rightness fluctuations as synchrotron cooling of high energy 



electrons (jUchiyama et al.ll2007l ). The 1-D probability distribution of most model parameters 
do not converge very well with multiple peaks except for those for a p , and E^, which are 
constrained by the 7-ray data directly, independent of other parameters in the hadronic 
scenario. 

Due to the lack of constraint from the 7-ray data, the parameters related to electron 
emission can not be well determined with relatively large la errors shown in Table [U, and 
there are strong correlations in the 2-D confidence contours, as shown in the middle panel 
of Fig. [2J The correlations among W e , E e c , and B are due to the facts that the synchrotron 
emissivity e oc W e B 2 E^ 2 , and the high energy cutoff of synchrotron emission v c oc BE e2 . v c 
is well constrained by X-ray observations, which leads to the anti-correlation between B and 
E e c . The anti-correlation between W e and B can be attributed to the measured luminosity 
of the synchrotron emission. The correlation between W e and E e c results from these two 
anti-correlations. The strong anti-correlation between B and a e is due to the radio to X-ray 
spectral shape, which, in combination with the anti-correlations between B and E%, and 
B and W e , leads to the correlations between a e and E e cl and a e and W e , respectively. The 
correlations related to 5 e can be attributed to the hard X-ray spectrum. The weak correlation 
between a p and E v c is due to the well-measured high-energy cutoff of the TeV emission. The 
strong anti-correlation between nisM and W p is due to the fact that the product of n ISM and 
W p determines the hadronic component of 7-ray emission. Therefore we can get a lower 
limit of W p according to the upper limit of nisM- The 2a upper limit of nisM is 0.01 cm -3 



lr The IC cooling is negligible compared with the synchrotron cooling for magnetic fields greater than 
10/iG. 

2 The break energy in the leptonic (and the following hybrid) model is very high due to the weak magnetic 
field, and the electron spectrum is identical to a single power law. 
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for T e = 10 7 K, corresponding to a lower limit of 1.6 x 10 52 erg for W p . The dependence of 
these results on T e will be discussed in subsection 2.4. The constraint on W p requires a total 
energy of CR protons much higher than typical energy outp ut of a supernova explosion, say 
10 51 erg ferlison et al.lboiol ; Izirakashvili fc Aharonianlboioh . 
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Fig. 2. — Same as Fig. [T] but for the hadronic scenario. The dashed line indicates the 
thermal emission for T e = 10 7 K and tiism = 0.011 cm" 3 . 

Table [2] shows clearly that the fits to GeV and TeV 7-ray data are significantly improved 
in the hadronic model. The average residual for the HESS data becomes ~ 1.14a - . The 
fit to the X-ray data also improves a bit due to a larger value of 5 e . Note, however, the 
reduced y 2 of the global fit is about 1.46, which is still too large to be attributed to pure 
statistical errors. There are some systematical effects in either the data or the model. In 
particular, contributions to the y 2 are dominated by the X-ray data in all these models, as 
shown in Table 2. It suggests that our simple one-zone synchrotron emission model does not 
give a sufficient description of the spatially integrated emission spectrum. Compared with 
the leptonic model, most reduction in the value of the \ 2 °f the hadronic model comes from 
improved fit to the 7-ray data. 

In the hadronic model the magnetic field is large, which suppresses the IC contribution to 
the 7-rays from energetic electrons and makes hadronic contribution to the 7-ray dominant. 
The strong magnetic field of the best fit model implies very efficient energy loss near the 
cutoff energy of electrons, which not only introduces a spectral break in the ov erall electro n 
distribution but can also render the high energy cutoff sharper with 5 e = 2.1 (|Blasill2010l ). 
It is interesting to note that the spectral index for electrons a e is consistent with that for 
protons dp, which is expected if the acceleration of these high-energy particles is due to the 
same physical process. The difference in their high-energy cutoffs can be attributed to the 
difference in the energy loss rate of protons and electrons near the cutoff energy. Therefore 
E e c should not be compared to E v e directly. 

To reduce uncertainties of the model parameters in the hadronic scenario, we also con- 
sider the hadronic model with the constraint that a e = a„. The 1-D probability distribution 



- 10 - 



and the 2-D confidence contours, and the best fit SED for such a model are shown in Fig. 
|3j Compared with Fig. [21 the probability distribution of the model parameters are better 
converged and the correlations in the 2-D confidence contours are weakened. The accept- 
able model parameter space is reduced significantly. The fitting results are also compiled 
in Tables [1] and [2] with a star mark. The value of the best fit model parameters agree with 
those of the hadronic model and the la error of parameters related the electron emission are 
reduced significantly. The values of the \ 2 are essentially the same as the hadronic model. 

The value of the spectral index for the best fit model is always much less than 2, a result 
difficult to accommodate with the diffusive shock model. This in combination with the low 
limit on W p poses one of the most serious challenges to the hadronic scenario in the context 
of diffusive shock acceleration of SNRs. We also note that the best fit electron spectral 
index a e is harder for the hadronic model than that for the leptonic model. Although the 
current radio data does not give a good constraint on synchrotron spectral index, a better 
measurement of the synchrotron spectrum in the future will be helpful in distinguishing these 
two scenarios for the TeV emission. 
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Fig. 3. — Same as Fig. |2]but for the hadronic scenario with the extra requirement of a e = a 



v 



2.3. Hybrid scenario 



The results above agree with previous studies (ITanaka et al.l I2008I ; lAharonian et al. 



20061 ). They demonstrate clearly the strengths and problems associated with the leptonic 
and hadronic scenarios. The leptonic model has fewer parameters, most of which are well 
constrained with the MCMC method by fitting the SED. The fact that it can give reasonable 
good fits to the overall SED with reasonable values of the parameters may be considered as 
evidence for such a scenario (jFan et al.ll2010al ). On the other hand, the relatively high values 
of x 2 , especially for the 7-ray data, suggest that the model may not be complete. It has 
been shown that the TeV emission fr om SNR RX J1713.7-3946 may have significant contribu - 
tions from energetic protons as well (jZirakashvili &: Aharonianll2010l : iKatz fc Waxmanll2008l ). 
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However, in the most general case, at least three more parameters need to be introduced 
to characterize the distribution of accelerated protons, which leads to strong degeneracy of 
the model parameter space. And the challenges to the hadronic scenario, namely hard spec- 
tra of accelerated particles and the lower limit on W p , does not appear to depend on this 
degeneracy in the regime of parameter space explored above. 

In the subsection 2.2, we demonstrate that the acceptable model parameter space may be 
reduced significantly by considering some physically motivated constraint on the accelerated 
particle distributions. We have argued that the spectral indices of accelerated electrons and 
protons should be comparable in the relativistic energy regime where the energy loss can be 
ignored. In the hadronic scenario, due to the presence of strong magnetic field, there is a 
spectral break in the electron distribution and the high energy cutoffs of electrons and protons 
do not need to be the same. However, in the leptonic scenario, the magnetic field is so weak 
that the energy loss does not affect the distribution of both electrons and protons. For these 
high energy relativistic electrons and protons, their gyro-radius only depends their energy 
and the magnetic field. We would expect that mechanisms of charged particle acceleration 
will lead to identical particle distributions exc ept their normalization , which is determined by 
different injection processes at low energies (IPetrosian fc Liull2004l ). To reduce the number 



of model parameters, one may therefore consider the hybrid scenario where a e = a p , 5 e = <5 



and E e c = E v c . As we will show below, this leads to a hybrid explanation to the high energy 
7-ray data. The total number of free parameters in the hybrid model is now only 7. 





Fig. 4. — Same as Fig. [T]but for the hybrid scenario. The density for the thermal emission 
indicated by the dashed line is 0.009 cm -3 . 

The model parameter distributions and best-fit SED are shown in Fig. HI The back- 
ground electron temperature is still adopted as 10 K. The best fit results and x 2 values are 
also listed in Tables Q] and [2j We see that the parameters of the electron component do not 
change significantly compared with the leptonic scenario primarily due to the dominance 
of TeV emission by relativistic electrons through the IC process. The hadronic component 
dominates the GeV 7-ray emission. Note that since E% = E%, the 7-ray spectrum of the 
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hadronic component cuts off at a lower energy than the leptonic component. Neutral pion 
decays into two 7-ray photos. The cutoff energy of the 7-ray spectrum is at least a factor 
of 2 lower than the cutoff energy of the corresponding proton distribution. For the IC emis- 
sion, the cutoff energy of the 7-ray spectrum can be the same as the corresponding electron 
distribution. The fit to TeV data shows improvement compared with the leptonic model, 
with average residual changing from ~ 2.3a to about 2.0a. However, the hadronic compo- 
nent seems to overproduce the GeV flux, resulting in an even larger Xcev than the leptonic 
scenario. The \ 2 value for X-ray data does not change significantly. The correlation of the 
model parameters are also similar to the leptonic model. And the strong anti-correlation 
between nisM and W p is still due to the fact that the observed emission is determined by 
the product of the two. Although the model has a weak magnetic field and relatively soft 
distributions of accelerated particles, the 2a lower limit of the proton energy of 1.0 x 10 52 
ergs is comparable to those of the hadronic models, which still challenges the energetics of 
the SNR. 



2.4. Dependence of T e 

The energy content of relativistic protons is poorly determined due to the high uncer- 
tainty in tt-ism- Perhaps the only observati on one can use to constrai n tiism is the lack of 



thermal X-ray emission from the remnant ( jCassam-Chenai et al.l 120041 ) . To derive a robust 



constraint on n^M, one however needs to consider the heating of electrons and the ionization 
of ions in the background plasma, both of which are not well understo od though a prelim - 



inary attempt has been taken to model these processes quantitatively ( lEUison et al.ll2010l ). 
Here we assume that electrons have reached ionization equilibrium with the ions and so that 
the Raymond-Smith code can be used to calculate the thermal emission. 

The results above do not differ significantly for different values of T e except for the 
constraint on n^M and accordingly W p . In Fig. [5] we show the results for the hybrid scenario 
with T e = 10 6 and 10 8 K. For T e = 10 6 K most of the line emission has energies lower 
than 0.5 keV, which is below the lower limit of the Suzaku data. However, the emission 
in the X-ray band is sufficient to lead to well-constrained uism ~ 0.2 cm -3 and W p ~ 10 51 
ergs. The model also predicts strong emission below the X-ray range. A significant thermal 
component also helps to slightly improve the fit to the X-ray data. The 2a upper limit of 
n^M for T e = 10 6 K is 0.2 cm -3 , which is much looser than that for T e = 10 7 K. For T e = 10 8 
K the 2a upper limit of n^M is 0.02 cm -3 , which is also higher than 0.009 cm -3 for T e = 10 7 
K. The 2a lower limits of W p are 3.5 x 10 50 , 1.0 x 10 52 , and 4.4 x 10 51 ergs for T e = 10 6 , 10 7 , 
and 10 8 K respectively. 
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Fig. 5. — Same as Fig. H]but for two other background electron temperatures: 10 6 (upper) 
and 10 8 (lower) K. The corresponding density nisM for the thermal emission is 0.2 cm -3 
(upper) and 0.02 cm -3 (lower). 




T„(K) T„(K) 



Fig. 6. — The 2a upper limits of tiism (left) and lower limits of the energy content of 
non-thermal protons W p (right) as functions of the temperature of background electrons 
T e for different models. The shaded region can be excluded due to heating of electrons by 
background protons through Coulomb collisions. 

To demonstrate how the constraint on n^M and W p vary with T e , we repeat the MCMC 
calculation for a series of T e . The 2a upper limit of nsiM and lower limit of W p are shown 
in Fig. |6l In the left panel the constraints on tt-ism of the three scenarios are shown, while 
in the right panel the constraints on W p are relevant for the hybrid and hadronic models. 
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The constraints on tiism are almost the same for the leptonic and hybrid models since their 
parameters for the synchrotron emission are similar. For the hadronic model the result is 
slightly different primarily due to difference in the electron distribution. The lower limit on 
W p for the hybrid scenario is a bit smaller than the hadronic model because in the hybrid 
model the IC component from the lepton population has a significant contribution to the 
7-ray. The shaded region in the left panel is excluded by considering the heating of electrons 
by ions through Coulomb collisions with the scaling relation of Eq. (1). In the right panel 
the shaded region is derived by applying the relation W p /10 50 erg ~ 1.3 cm~ 3 /?iisM! which 
is the approximate relation required to reproduce the high energy 7-ray emission for the 
hybrid and hadronic scenarios. It can be seen that in general W p needs to be greater than 
~ 4 x 10 51 erg for the hybrid model. For the hadronic scenario the requirement of W p is 
even l arger. Such a large value of proton energy seems to be unacceptably high for a typical 
SNR ( lEUison et al.ll2010t IZirakashvili fc Aharonianll2010f ). The 2a upper limit o f tt-tsm is less 
than 0.03 cm -3 , which is consistent with that obtained by lCassam-Chena'i et all (120041 ) from 
XMM-Newton observations. 



Discussion and Conclusion 



Since the discov ery of synchrotron X-ray emission from the forward shock of SN 1006 
( IKoyama et al.lll995l ). it has been established that shocks of SNRs can accelerate electrons 
to tens of TeV. The detection of TeV emission directly from a few s hell type SNRs confirmed 
this conclusion ( iLazendic et al.ll2004l ; lAharonian et al.ll2005l 120081 ). However, the nature of 
the TeV emission is still a matter of debate. X-ray observations also reveal strong thermal 
emission in the interior of a few young SNRs. Thermal X-ray emission from the forward 
shock, on th e other hand, is u sually weak and only detected in relatively old remnants, such 
as RCW 86 (IVink et al.l 120061 ) . Since the observed TeV emission is well correlated with the 
forward shock, the lack of thermal X-ray emission implies low gas density and therefore high 
energy content of accelerated protons in the hadronic model for the TeV emission, which is 
difficult to accommodate with the SNR theories. The leptonic model in general gives a poorer 
fit to the TeV spectrum than the hadronic model. More detailed modeling is necessary to 
distinguish these emission models. 

Using the MCMC method, we systematically investigate the parameter space of models 
for the multi-wavelength emission of SNR RX J1713. 7-3946. The high quality of observa- 
tional data, especially X-ray data from Suzaku and TeV 7-ray data from HESS, enables us 
to get very good constraints on most model parameters and better understand the emission 
mechanisms. The radio and X-ray emissions are thought to be produced by the synchrotron 
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radiation of relativistic electrons accelerated in the SNR. The high energy 7-ray emission 
(from GeV to TeV) can be produced through the IC radiation of electrons scattering off 
background low energy photons, and/or the decay of ir° generated through CR proton-ISM 
collisions. We study three kinds of scenarios, the leptonic one, hadronic one and a hybrid one, 
distinguished through the emission mechanism of high energy 7-rays. Thermal emissions, 
including continuous and line emissions, are included in the modeling to constrain density 
of the background plasma with the absence of thermal emission in X-ray observations. 

The global fit of these three scenarios shows that: 1) the goodness-of-fit is the worst 
for the leptonic model and the best for the hadronic model; 2) the X-ray data can set an 
upper limit on the background ISM density, which is 0.03 — 0.009 cm -3 depending on the 
temperature of the background electrons; 3) the upper limit of tiism leads to a lower limit of 
energy content of relativistic protons W p > 4 x 10 51 (6 x 10 51 ) erg for the hybrid (hadronic) 
model, which seems to be too high for a typical core-collapse supernova. The well constrained 
parameters of the leptonic scenario are more physically acceptable, although the hybrid and 
hadronic scenarios give smaller x 2 values. In some physic ally motivated m odels of electron 



acceleration in SNR, the goodness-of-fit can be improved (|Fan et al.ll2010bl ). and our overall 



results appear to favor the leptonic scenario (lEllison et al.ll2010l ). However, systematic errors 



in the X-ray and TeV band are significant in all the emission models studied so far, which 
demands more advanced modeling. For example, the observed source structure is much 
richer than the simple uniform zone assumed in all the emission models. On the other hand, 
given difficulties in X-ray and TeV observations, the errors of the relevant data may also be 
underestimated, especially the HESS fluxes at the low energy range, which doe s not appear 



to sm oothly match the recently obtained GeV fluxes with the Fermi observatory (lAbdo et al. 



20111 ). 



In the following we discuss several possible alternatives of the hadronic scenario, in 
order to provide an independent judge of the price needed to pay to keep the hadronic model 
working. The most serious challenge to the hadronic and hybrid models is the high energy 
content of relativistic protons inferred from the upper limit of the background density. This 
high energy content is in excess of typical supernova explosions, but could be explained by 
hypernovae. The total kinetic energy of the SNR with a shock speed U may be estimated as 

27T „ 3 __ ,. .. ,„.-,,) ( R \ 3 ( U \ 2 f n ISM 



K = y «W - 6 x 10- j j {^^) ergs. 

With the current best estimate of the distance at lkpc, the radius of the remnant is 10 pc. 
For such a remnant size, the total kinetic energy content (K) of the remnant is much less 
than typical hypernovae and less than the energy content of relativistic protons. However, 
the distance to the SNR is not well determined. The lack of thermal X-ray emission di- 
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rectly constrains the integration along the line-of-sight of the thermal emissivity, which is 
proportional to n 2 SM R. The radius R and shock speed U of the SNR scales linearly with the 
distance D. The upper limit on nisM therefore is proportional to D~ x l 2 \ K is proportional 
to D 9 / 2 . The 7-ray flux produced via hadronic process is proportional to W p nisM/D 2 . The 
lower limit on W p therefore is proportional to D 5 ^ 2 . The total kinetic energy of the SNR 
decreases more rapidly than W p , making the discrepancy between these two energies worse 
for shorter distance D. An increase of D by a factor of 3, on the other hand, can make the 
upper lim it on W p less than K , which is now ~ 8 x 10 52 erg, at the high end of observed hy- 
pernovae ( iNomoto et al.ll2005l ). This may explain the rareness of this kind of shell type TeV 
SNRs with non-detection of thermal X-ray emission. For the fe w tens of SN Rs observed in 
X-rays, two of them show complete absence of thermal emission ( )Vinkll2006l ). The efficiency 



of pro ton acceleration also needs to be greater than 50% for this model to work ( iHelder et al. 
20091 ). The upper limit on the density is then about 0.02 cm -3 , which is not too different 
from values obtained above and from previous studies. 

Alternatively, the hadronic and hybrid models can work by overcoming the upper limit 
on the ISM density placed by the upper limit on the thermal electron emiss i on. T his limit 
arises from using the standard ion-electron coupling term from lHughes et al.l ( 120001 ) produc- 
ing a lower limit on the electron temperature (Eq. (1)). But this assumed a standard profile 
for the remnant. If the emission of this remnant is instead produced by the interaction of 
the supernova shock striking a thin shell of material ejected from the star prior to collapse 
(e.g. an outburst from a luminous blue variable or binary mass ejection), the electrons might 
still be cold. However, to do so would require very fine-tuned arguments (for the electron 
temperature to be below 2 x 10 6 K, the shock must have hit the stellar ejecta less than 25 yr 
ago). Even with this fine-tuning, this thin shell of material would not have a large enough 
emitting volume to explain the observed high-energy emission. A thin shell explanation 
seems unlikely. It was also suggested that the effective particle acceleration would make 
the post-shock medium be less heated (IDrury et al.l 120091 ). however, the extreme condition 
in such a scenario seems difficult to satisfy, and the no n-linear effects will also make the 
medium be heated more significantly (lEllison et al.ll2010l ). 



A plausible explanation of the low inferred density of the background plasma is the 
assumption of ionization equilibrium adopted in calculating the thermal emission, which 
however may not be valid. For the case with a large fraction of neutral medium, n e <C uism, 
the constraint on n^M from the thermal emission, which is proportional to n e riisM, will be 
much weaker than the case with fully ionized gas. Therefore the requirement of W p can be 
smaller and the difference between W p and the power of typical supernova can be reduced. 
In such a case, the low density of ionized ISM implies that neutral gases can penetrate a 
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depth of 



6 x 10 
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u 



0.03 cm-V V 4500km / S y V10 7 Ky 



(■ 



0.23 



cm 



before being ionized by the free electrons ( iChevalier fe Raymondlll978l ). For very low density 
of ionized plasma in the shocked downstream region, this depth can be a significant fraction of 
the radius of the r emnant, and one may expect strong H Q emissions (jGhavamian et al.ll2007l ; 
Helder et al.ll2009l ). Better treatment of the ionization balance in the shock downstream may 
address this issue. 

Finally we mention that neutrino signal can be used to test the hybrid/hadronic model 
of the 7-ray emission. It is shown that if the TeV 7-rays are predominantly produced by 
hadronic interactions, the accompanied neut rino signal might be det e cted in the up-coming 



km 3 neutrino detector, such as KM3N ET (IKistler fe Beacoml 120061 ; lYamazaki et al.1 12009 
Morlino et al1l2009bl ; lYuan et alJboioh . 
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